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Abstract. Simplified forms of the space-time discontinuous Galerkin (DG) and 
discontinuous Galerkin least-squares (DGLS) finite element method are developed 
and analyzed. The new formulations exploit simplifying properties of entropy en- 
dowed conservation law systems while retaining the favorable energy properties 
associated with symmetric variable formulations. 


Introduction 

The high-order accurate numerical solution of systems of nonlinear conserva- 
tion laws remains a computationally expensive endeavor. This article consid- 
ers simplified forms of the Discontinuous Galerkin (DG) and Discontinuous 
Galerkin least-squares (DGLS) finite element methods tailored to systems of 
first-order nonlinear conservation law’s with convex entropy extension. Cen- 
tral to the development is the Eigenvalue Scaling Theorem which character- 
izes right symmetrizers of entropy endowed systems of conservation law’s in 
terms of scaled eigenvectors of the corresponding flux Jacobian matrices. This 
yields a simplification of the DG and DGLS methods without sacrificing the 
pleasing energy properties of symmetric variable formulations. The next sec- 
tion briefly review’s a number of results in symmetrization theory discussed 
in detail in Barth [2.1]. 

1 Brief Review of Symmetrization and the Eigenvector 
Scaling Theorem 

Consider a system of m coupled first-order differential equations in d space co- 
ordinates and time which represents a conservation law process. Let u(x,£) : 
JR d x 1R + i-4 lR m denote the dependent solution variables and f(u) : IR m 
IRmxd t ^ e ^ ux vector ^e prototype conservation system is then given by 

+ = o (i.i) 

with implied summation on the index i. Additionally, the system is assumed 
to possess an scalar entropy extension. Let U{u ) : IR m ^ 1R denote an 
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entropy function and F(u) : IR m IR^ the entropy flux such that in addition 
to (1.1) the following inequality holds 


u.t + F* x . < 0 

(1.2) 

with equality for smooth solutions. In symmetrization theory for first-order 
conservation laws, one seeks a mapping u(v) : IR m IR m applied to (1.1) 

so that when transformed 

U.vV.t + f[ v V ]Xi = 0 

(1.3) 

the matrix u_ v is symmetric positive definite (SPD) and the matrices /* are 
symmetric. Clearly, if functions U(v),J ri (v) : lR m ►-> IR can be found so that 

n T =u v , ( f‘) T = r v 

(1.4) 

then the matrices 


u v = f v — J~] v v 

(1.5) 

are symmetric. Since v is not yet known, little progress 
Introducing the following duality relationships 

has been made. 

U(u) = v T (u)u - U(v(u )) 

(1.6) 

F‘(u) = v T {u) /*(«) - ^(u(u)) 

(1-7) 

followed by differentiation yields 


L'.u = v T 4* u T v yU - U yV v yU = v T 

(1.8) 

Cu = V T f' u + (/OV- - F v V, u = V T f‘ u 

(1.9) 


Equation (1.8) gives an explicit expression for the entropy variables v in 
terms of derivatives of the entropy function U(u) 

V T = [/« - (1.10) 


Finally, we require convexity of U(u) which insures positive definiteness of 
«,u and u,„ and implies hyperbolicity of (1.1) [5,9], viz., that the linear 
combination / u (n) = n, f' u has real eigenvalues and a complete set of real- 
valued eigenvectors for all nonzero n € IR d . This result follows immediately 
from the identity 

(u. v )~ l/2 f U (n)(u, v ) l/ - = («,„) _1/2 / „(n)(u.„)~ I/2 

^ y 

symm 

which shows that f u (n) is similar to a symmetric matrix. 
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1.1 The Eigenvector Scaling Theorem 

Next, we consider an important algebraic property of right symmetrizabie 
svstems which is used later in the implementation of the DG and DGLS 
schemes. Simplifying upon the previous notation, let Aq - u v * = f 
A { — AiA 0 and rewrite (1.3) 

+ A t V' £i — 0 . ( 1 - 11 ) 

The following theorem states a property of the symmetric matrix A t sym- 
metrized via the symmetric positive definite matrix To* 

Theorem 1.1 (Eigenvector Scaling). Let T G !R rlxn be an arbitrary di- 
agonalizable matrix and S the set of all right symmetrizers: 

S = {B G lR nxrl | B SPD, AB symmetric}. 

Farther . let R £ IR n * ri denote the right eigenvector matrix which diagonalizes 
A 

.4 = RAR~ l 

with r distinct eigenvalues, .1 = Diag(Ai I mi x mi • x m 2 * ■ • • .A r / mr . x m r )* 

Then for each B € S there exists a symmetric block diagonal matrix T = 

Diag r mf xmj that block scales columns of R, R = 

RT. such that 

B = RR t . A = RAR~ l 


which imply that 


AB = RAR t . 


Proof. Omitted, see [2]. ^ 

Note that this last formula states a congruence relationship since R is not 
generally orthonormal and A does not represent the eigenvalues of AB. The 
Eigenvalue Scaling Theorem is a variant of the well known theory developed 
for the commuting matrix equation AX — XA = 0, A,X € IRT*", see for 
example Gantmacher [6]. Examples of the Eigenvector Scaling Theorem for 
the Euler and magnetohydrodynamic equations are given in [2]. 


1.2 Generalized Matrix Functions with Respect to the 
Riemannian Matrix Ao 

For use in later developments, it is useful to define the following generic 
matrix function f(A) with respect to the Riemannian matrix A 0 

f Ao (A) = A 0 f(A^A). (1.12) 

This definition reflects the following steps: (1) multiplication of the system 
(1.11) by in order to restore a Euclidean metric, (2) invocation of the 
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matrix function on the matrix product Aq 1 Ai, (3) multiplication of the result 
by _4o to restore the original metric matrix. Proposition 1.2 shows that this 
generalized matrix function is symmetric and has a rather simple construction 
for symmetrizable systems by exploiting the Eigenvalue Scaling Theorem. 

Proposition 1.2. Barth [2Aj. Let T () denote the SPD right symmetn.zer of 
A such that A = .4.-1-,. .4 0 = RR T . and .4 — RAR ~ l . The generalized matrix 
function f ^ (A) is symmetric and defined canonically in terms of entropy 
scaled eigenvectors as 

f- Xo (A) = Rf(A)R T . (1.13) 


In later sections, the generalized matrix absolute value function |.4| ^ will be 
required 

1-41.4,, = m\R T ■ (i.i4) 

This matrix absolute value function has a natural generalization to IR d using 
an Lp-like norm definition 


• 4 1 





(1.15) 


which has a particularly simple form when p = 1 which is used later in the 
least-squares term appearing in the DGLS method 


d 

-4|i..4 0 = ■ (116) 

i= L 


2 Simplified DG, DGLS, and GLS Finite Element 
Methods 

Let i? denote a spatial domain composed of nonoverlapping elements T*, 
J? = UTi, T t nTj = 0, i ^ j and I n =]t n ,£ n+1 [ the n-th time interval. 
It is useful to also define the element set T = {7\, • • • , T\t\} an ^ edge 

set £ = {ei.e 2 ,.. . . e fj}. To simplify the exposition, consider a single vari- 
ational formulation with weakly enforced boundary conditions. By choosing 
the correct space of functions (discontinuous or continuous) and omitting the 
least-squares variational term, one can switch from the DGLS formulation to 
the DG or GLS formulations. In the GLS formulation [7,10], functions are 
continuous in space and discontinuous in time 

v* = {v ft |r h 6 (c°(/2xn) m , v f rx/ „ e (p*(Tx /"))”*} 

where v denotes the entropy variables for the system. In the DG and DGLS 
formulations [8.3], functions are discontinuous in space and time, i.e. 

v* = 6 (p*(Tx /"))“} . 
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Consider the prototype hyperbolic system for the space-time domain ftx [0, T] 
with boundary data g imposed on F via admissibility condition 

u t + f l Xt = 0 in ft 

*4“ (to) (g - u) — 0 on r (2.1) 

or in symmetric quasi-linear form for smooth solutions 

AqVj 4- = 0 in ft 

(n) (g - v) ~ 0 on r (2.2) 

with .4 (ti) = n, A, and A(n) — n t .4,. The combined GLS and DG schemes 
are defined by the following stabilized variational formulation: 

Find v h £ Y h such that for all w h € V* 

B(v h .w h ) OAh + B(v h ,w h ) L s + B(v h ,w h ) BC =0 (2.3) 


B(v.w)c,al=J J (-!*(») W.t - f l (v) -w, Xi )dxdt 

+ [ (ti>(f" +l ) • -u(v(<" +1 )) - w(t+) «(w(t"))) dx 

J Q 

+ [ ^ f(w(x-) - w(x+)) ■ h{v{x-),v(x+):n)dxdt 
B{v. W LS = + A t w Xi ^ • r ^.4qV.£ 4- A t v 


dx dt 


B(v,w) bc = J J w ■ h(v,g;n) dx dt 

where h denotes a numerical flux function and r a small m x m SPD matrix 
for the least-squares term. For theoretical and practical reasons, two numer- 
ical flux functions are considered. Both are of the form 

h(v-,v+-,n) = 1 ( f(v-;n ) 4- f{v+;n)) - ^h d (v-,v+;n) (2.4) 

and consistent with the true flux in the sense that h(v,v;n) = f(v;n). 


1. Symmetric Mean- Value Flux. This flux is motivated from the non- 
linear energy theory of Sect. 2.2. Define the parameterization v(0) = 
v(x_) + 0 [v]l + . The symmetric mean- value flux is then given by 

hg MV (t>_,w + ;n) = |i(u_,u + ;n)| S MV [v} x x * 

with 

|.4(w_,r+rn)|sMV = [ |.4(u(0); ra)|^ o d.8 . (2.5) 

Jo 

By construction, the matrix \A(v- , v+ : to)|smv is symmetric positive 
semi-definite. Using this form of flux dissipation (2.5), nonlinear entropy 
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norm stability of the DG, DGLS, and GLS formulations can be shown as 
discussed in Sect. 2.2. In addition, let 

Mn ) smv = f A(v(0):n)d0 ( 2 . 6 ) 

Jo 

denote the integral mean-value matrix for which the following useful prop- 
erty exists 

[/(«)]* t = -4(n)sMv ( 2 - 7 ) 

which is a necessary ingredient for optimal discontinuity resolution. To 
prove stability of other (more practical) forms of flux dissipation, one 
formally needs only show that the new form is more energy dissipative 
than the symmetric mean- value form in the following sense: 

[v]i:-hi uv <[v}i:.h d . 

2. Discrete Symmetric Mean- Value Flux. The discrete symmetric mean- 
value flux function replaces the state-space path integration in (2.5) by 
.V point numerical quadrature 

^DSM \{v-,v+-,n) = |.4(u_,u + ;n)| D SMV [ v ]*+ 

with 

jV 

|.4(v_. n)|osMV M* t = ;c«^w,);n)u o k: ( 2.8) 

q=i 

where w q and 0 q denotes the quadrature weights and positions. In form- 
ing this flux, recall from the Eigenvalue Scaling Theorem that |A|^ o = 
R\A\ R t . This flux function is of practical interest since it is easily formed 
and has a relatively straightforward Jacobian linearization as will be 
shown later. The absolute value in this equation renders the state-space 
path integration from t/_ to v+ slope discontinuous whenever entries in A 
change sign. In this case, to retain accuracy of the numerical quadrature 
at n < m possible points of slope discontinuity, the path integration is 
further subdivided into subintervals, e.g. [v_,v+] = [v_, v\] U [vJjVj] U 
. . . U [v n ,v+] where v * = v(0i) is a location 9i such that an entry of 
.1 vanishes. In practice, satisfactory results [1] have been obtained using 
2-point Gaussian quadrature rules (which integrate cubic polynomials 
exactly). 


2.1 Linear Energy Analysis 

Due to length constraints of this article, we simply restate some relevant 
theorems given in Barth [2,1] concerning energy boundedness of variational 
form (2.3) for systems of hyperbolic equations. 
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Theorem 2.1. Global Energy Stability (Linear Hyperbolic System), 

The variational formulation (2.3) for linear hyperbolic systems is energy sta- 
ble (modulo data g j with the following global energy balance: 


i^,-, Q ~ 2 :l- 4 ou./ + ■4i».x 1 |lr.flx/- +H([ U ]x-):.4 , f x/- ) 

: 0 \ / 

V — l /V-l 

£ \ v ( t - Jlldo.r? = )ll.4 0 ,r? + £ 2 ( v '5){-d- ).r x/» 


a=0 
V-l 


n=0 


n=0 


Proof. Omitted, see 2.1]. E 

This energy balance equation formally bounds the energy at time in terms 
of initial data and inflow boundary data. Next, we consider the full nonlinear 
system of conservation laws. 


2.2 Nonlinear Energy Analysis 

Before presenting the nonlinear energy result, we prove a general lemma and 
consequential corollaries concerning entropy function/flux jump identities at 
space-time slab interfaces. Note that throughout this section, we utilize the 
state-space parameterization 

v(0) = v(x- ) + 0 [v] x x + 


(similarly across time slab interfaces) for use in state-space path integrations 
and the interface averaging operator 

Lemma 2.2. Interface Jump Identities. Let Z(u),Z(i>) : IR m 1R be 

twice differentiable functions of their argument satisfying the duality relation- 
ship 

Z(u) + Z(v) = Z, v v . (2.9) 


The following jump identities hold across interfaces 

[Z]^ - [Z. V ] X S v(x+) + f\l - 0)[v]l± Z^ v {vm Wit d9 = 0 (2.10a) 

Jo 

[Z}*+ - [Z, v ]* x + *(*_) - [ 1 0 [v] x x +- Z v , v ( v(9)) [v]lt d0 = 0 . (2.10b) 

Jo 

Proof. Recall the following forms of Taylor series with integral remainder 

[Z]'- - Z.„(x.) [t>]' + + I' 9 Z. v . v (v(9)) [v]*+ d0 = 0 (2.11a) 

Jo 

[Z]' + -Z,|I.) [v] x x - - f\l-0) [v} x x +-Z' V , v (v(0)) Wx! d« = 0 ( 2 . 11 b) 

Jo 
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and the jump form of (2.9) 

[-ix! - ( Mx! + «»»;: . (2.12) 

Combining 2.11a). 2.11b) and (2.12) yields 

[z} J t :-((v) i x :[ 2 . v r i :+^J\i- 29 ) {v\^-z. v . v {vm «: m = 0 . ( 2 . 13 ) 

Finally, algebraically manipulating this form together with the mean-value 
identity 


• 7 l x + 


/ 


Z.v,v(v(9)) d0 


produces the stated lemma. 


(2.14) 

□ 


Corollary 2.3. Temporal Space-Time Slab Interface Identity. Let t ± 

denote a temporal space- time slab interface. The following entropy function 
jump identity holds across time slab interfaces 

f n (WtZ -v r (4)M«!) dx + ^111 Ill 2 4 0 .n = ° (2-15) 

where 

ill Mtl 1, n= f [ 2(1 — 0) [v]j + • .40(17(0)) [i>]J + dddx > 0 . (2.16) 

J qJ o 

Proof. Set Z — U . Z —U with U, v = u T and v ,v = *4o in Lemma 2.2 and 
replace with t± in (2.10a), see [10] for an alternative form. □ 


Corollary 2.4. Spatial Space-Time Slab Interface Identity. Let x± 
denote a spatial element interface. The following entropy flux jump identity 
holds across spatial element interfaces 


d9 = o . (2,i7) 


Proof. Set Z = F'. Z = i = 1, . . . ,d with T' v = ( f i ) T and T' v v = .4, 
in (2.13) of Lemma 2.2. □ 


Note that in actual numerical calculations, it is desirable to use the variational 
form given by (2.3) since integration by parts has been used to insure exact 
discrete conservation even with inexact numerical quadrature of the various 
integrals. For analysis purposes, however, it is desirable to use the following 
equivalent non- integrated-by= parts formulation: 

Find v h E V h such that for all w h E V h 

B(v h .w h ) gal + B(v h . w h ) ls 4- B(v h ,w h ) B c = 0 


(2.18) 



Simplified DG Methods for Systems of Conservation Laws . . . 


9 


Biv. ) G al = / [ w ■ (u, + f l x .{v)) dxdt 
J i"Jn 

4- J u>(t" ) • [ujjj dx 

+ j ^ J' ^ ML ' h' l {v{x-),v(x+)\n)dxdt 

+ J t ' [ /( u:n )lL 

S(v. w.Ils = /X/ t m t 4- • r (^AoV't 4 dxdi 

B(v.w)bc = J J w * \ (/(^;n) - f{v:n) - h d {v.g:n)^ dxdt 

where h d denotes the flux dissipation term incorporated into the total nu- 
merical flux. 

Theorem 2.5. Global Entropy Norm Stability (Nonlinear Hyper- 
bolic System). The variational formulation (2.3) for nonlinear systems of 
conservation laws with convex entropy extension and symmetric mean-value 
flux dissipation 

hsMvfv-. v+:n) = |-4 |smv » I^Ismv = [ \Mv{0)\n)\^d0 

Jo 

is entropy norm stable (modulo data g) with the following global balance: 

N-i / 

( IIIMln ||| 2 4 0 ^ ■+■ 2\\AoV,t + AiV.xM r,!7x/" + 5^(M^t)jU| ie x/ 

ri=0 V 

rx/ .+ [2U(t?)dx= f2U(t°_)dxfjj (v,g) { _ r)r ^+G n r (g.v)) 

n=0 }Q J ° n= 0 

with 

|i(n)| = f 2(1 -9) (A + (v(6y,n) Ao -A-(v(l-6y,n) A ) d6 
Jo v 

and 

Gr(g.v) = J J (F(g;n)- J 9 g ■ A(v{6)\n) g do'j dxdt . 


Proof. Construct the energy balance for the interval [t® s t!^] = U ^_ 0 1 I n 
by setting w = v and evaluating the various integrals. Consider the time 
derivative integral 


/ I v T u t dtdx= I / U t dt dx 
JQJi n 




dx 



in 
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and combine with the jump integral across time slabs. From Corollary 2.3 

Inh vFu dt dX ~ in dx = J n dx + §HI W/t lll^ 0if? ■ 

When summed over all time slabs, the first term on the right-hand-side of 
this equation vanishes except for initial and final time slab contributions. 
Next, consider the spatial operator term and apply the divergence theorem 


[ [ v T f l Xt dxdt = f f F l x dxdt 

J i"Jn ' Ji*Jn 

■la — [F(v;n)]*^ dx dt + J J F{v\n)dxdt 


where F(v : n) = n. F l (v). From Corollary 2.4 and the definition of |*4|, it 
follows that 


f l J n vT f'*> dxdt ~ j t '[/( n )]x I + \ Wx! h SMv)dxdt 

= Y2 + j'J r F {v-,n)dxdt . 


In summary, collecting terms we have 


B(v, v)g.\l 


= f [U\r 

J n 


dx + -| 




— J J F(v:n)dxdt 


ee£ 


A\,exl n 


The least-squares integral produces a pure quadratic form without modifica- 
tion 


£(tM>) LS = idoV.t + Av, x ,\\t nxI n . 


Finally, consider the boundary condition terms and apply Corollary 2.4 


B(v.v)bc = J J - f(v;n)) - - |.4(n)|s\iv(ff - «)) dxdt 

~ f J r (B{g;n)-F(v,n)+ ^ J{l-20)(g-v)A(v{0)-,n ) (g-v)dd 

• (/(s; n) - f(v, n)) - t v • |.4(n)| S Mv(9 - t»)) dx dt 

~ f f n ) ~ F ( v > n ) ~ J e 9 ■ Mv(0)) 9 dO 

• |-4(n)| v ~ v -T (n)g^jdxdt . 


Combining the above results, summing over time slabs, and multiplication 
by two yields an entropy norm balance equation which bounds the global 
entropy norm of the system at the final time T in terms of the initial data 
and boundary data g. □ 
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Remark 2.6. Note that when the A t matrices are assumed constant, /* = 
A t v and F* = so that the additional term G r p{g.v\n) vanishes 

identically and a one-to-one correspondence of terms between theorem 2.1 
and theorem 2.5 is achieved. 


2.3 A Simplified DG Method in Symmetric Form 

DG Flux Formulas. Simplification of the discontinuous Galerkin method fol- 
lows by choosing the discrete symmetric mean-value flux function proposed 
earlier, i.e. 

1 V 

h D ^i\ (v-.v+:n) = I (/(u_;n) + f{v+; n))-- ^ u.',|.4(u(0,)); n)| 4o 

< 1 = 1 

with \A\^ = R\A R t . v{0) - v(x-) + B >' usin S sufficient order 

numerical quadrature and subdivision of the state-space path integration 
at points of non-differentiability, the hoswv ^ ux can be rna( ^ e arbitrarily 
close to h^ NIV for which nonlinear stability in the DG method follows from 
the analysis of Sect. 2.2. Suppose that elements of \A(v{O q )); n)| remain 
bounded for 6 £ [0. 1] independent of N . Using N point Gaussian quadrature 

'^DSMV ” ^SMvIb = 0([u] ) - 

Next, we consider single-point quadrature formulas. 

Theorem 2.7. Discrete Symmetric Mean- Value Flux. Letv m be a state 
such that 

[»]': • |.4(u. : n)| [v]* + = sup [v]*+ ' |.-i(»(0);n)lj o M*! ■ 

O<0<1 

The variational formulation (2.3) with numerical flux function 

h DS \i\.{v-,v+:n) = \ (/(«_;n) + f(v+:n)) - l|-4(w.); n)|^ o [v] x x + _ 

(2->9) 

is energy bounded in the sense of Theorem 2.5. 


Proof. It is sufficient to show that the given flux dissipation 

^dsmv* — |-4(u.;n)| /io [ v ] x t 

exceeds the symmetric mean- value value flux dissipation. This is reflected by 
the algebraic condition 

[ v ] xt * ^SMV < Mx* ' ^DSMV* ■ 

From the symmetric mean-value flux definition 

K: ■ /»smv = [v\iz • f o i-w);n)u o do [«]^: 

< sup HV-|.W);n)|, 0 [<: 

0 <#< 1 
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= [<M.4( V .;n)|^ o [v]* x + 

= Mxt ' ^DSMV ■ 

This establishes nonlinear stability of the DG method using the simplified 
Hux function. □ 

Remark 2.8. Unfortunately, the state v m is not generally known in closed 
form. Cockburn and Shu [4] have shown impressive results using the simpler 
Lax- Friedrichs flux. It is straightforward to derive a corresponding “symmet- 
ric Lax-Friedrichs" numerical flux function 

hsLF(v-. u + : n - ^ (/(*?_: n) 4- /(v+;n)) - ^A rnax A 0 (v.) [v]* + 

with A max = sup 0<ti<1 (maxi<j< m (T z , (v(0)))). Nonlinear entropy norm sta- 
bility follows starting from Theorem 2.7 

[«;*! h SM\- < Kxl ■ n)U o [»]{£ 

= ■ R{v,\n ) |.l(u„;n)| R T (v,;n) [ 1 ;]^+ 

< sup ( max (.lu(u(0)))) [t>]*+ ■ R(v.;n) R T (v,;n) [v] x x + 
0<*<1 M<Km 7 

= A max [v} x x + _ ■ .4 0 (u.;n ) [«]*+ . 

Finally, for systems such as the Euler equations of gas dynamics that exhibit 
the property ^^f~ z,ZjZ k zi > 0, |z| / 0, we have 

[ v Yi~_ --4o(u.: n) [w]^ < max([t>]'+ -.4 o (t?(0); n) [v\l + _ , [«]*+ -.4 0 (t;(l); n) [»]*+) 
thereby avoiding the need for calculating v * altogether, see [1] for details. 


DG Jacobian Derivatives. Using the discrete mean- value fluxes, it becomes 
straightforward to compute Jacobian derivatives of various terms. For exam- 
ple. to compute derivatives of |-4|^ with respect to a vector w, chain-rule 
differentiation is used 


d\A(n)\ Ao dR(n) 


A{n)\R J ( n)+R{n ) 


d\A(n)\ 


fl r (n)+fl(n)|.l(n)| 


dR T (n) 


dw dw dw V-/ dw 

Note that a high degree of computational efficiency can be achieved in the 
calculation of these Jacobian terms by exploiting the transpose symmetry of 
intermediate products. 


2.4 Simplified Least-Squares Stabilization in Symmetric Form 


Consider an isoparametric element mapping J hx from a unit element space 
£ to a physical space x. In the papers by Hughes and Mallet [7] and Shakib 
[10], they proposed the follov/ing form for r appearing in (2.3) on a mapped 
element 


/ d \ Up d 

= ISIplo- l*Uo = E \ B ‘\ P B ‘ = • (2-20) 


\i = o 


j = o 
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Equation r2.20) is of the same form given earlier in (1.15). In standard im- 
plementations of least-squares stabilization, p — 2 is used. In light of the 
Eigenvector Scaling Theorem 1.1. it is useful to revisit the derivation of r 
with p = 1. Let B l — B l To, from (1.16) it follows that 

d 1 _1 

|V*°|.4o + Y. I V CI R( n ‘) |-l(n')i R r n‘) 

t= 1 

using the entropy scaled eigenvectors R(n l ) of B l . This represents a substan- 
tial simplification of the r matrix calculation. 

3 Concluding Remarks 

Simplified forms of the DG. DGLS, and GLS schemes have been presented 
and analyzed for first-order systems of conservation laws with convex entropy 
extension. Numerical examples are given in [1] using linear, quadratic, and 
cubic element approximation. 


ri = \B\ 


i.A„ 
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